AO-A054  677  MISSION  RESEARCH  CORP  SANTA  BARBARA  CALIF  F/«  20/14 

PROPAOATION  OF  MICROWAVE  SATELLITE  COMMUNICATION  SI6NALS  THR0U«->ETC(U) 
SEP  77  R W HENDRICK  ONA001-76-C-013S 

UNCLASSIFIED  MRC>R-334  DNA>4412T  NL 


flU  NO.:::^ 

DC  FILE  COPY  AD  A 05 467  7 


r 


flHl  LOHiHER  TRAN 

DNA  4412T 

PROPAGATION  OF  MICROWAVE  SATELLITE 
COMMUNICATION  SIGNALS  THROUGH 
STRIATED  MEDIA 

Mission  Research  Corporation 
735  State  Street 
Santa  Barbara,  California  93101 

September  1977 

Topical  Report  for  Period  November  1975— September  1977 
CONTRACT  No.  DNA  001-76-C-0135 


THIS  WORK  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDTSeE  RMSS  CODE  B323076464  S990AXHB05405  H2590D. 


D C 


UNCLASSIFIED 


security  classification  of  this  page  (When  Dmte  Enleted) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


|2.  GOVT  ACCESSION  NO.  3_  RECIPIENT'S  CATALOG  NUMBER 


-|-|-| t ■■■  ■ 

"'propagation  of  JUCROWAVE  ^TELLITE pOMMUNICATIONyroplcal  /epmrt 

■>»T/-»»TAT  C*  'PTlT>rVTT/^ll  C»  O'Ti  T A 'ri7TV  ^T?T\TA  ^ H C.  . C 


&/>ERlOO  COVERED 


J^GNALS  THROUGH  glRIATED^EDIA- 


/■ 


7.  AuTHORr*; 


^ W. ^Hendrick,  Jr.  j 


(lE^ 


Nov  75“Sep  77, 

■*^Tr^r.-T  HUMBER 


"^OR  GRANT  NUMBERCs; 


DNA 


-76-C7<Ji35  ' 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Mission  Research  Corporation''' 


735  State  Street 

Santa  Barbara,  California 


y / 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  4 W^fOfvUNIT  NUMBERS 


NWED  Sub  tasty  S99QAXW054-05 


93101 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Director 

Defense  Nuclear  Agency 
Washington.  D.C.  20305 


<4.  monitoring  agency  name  a A0DRESS('K  diHerent  from  Controlting  OHice) 


15.  SECURITY  CLASS 

UNCLASSIFIED 


ISa.  DECLASSIFICATION  DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  ST  AT  EM  EN  T ('of  this  ReporO 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  from  Report) 


18.  supplementary  notes 


This  work  sponsored  by  the  Defense  Nuclear  Agency  under  RDT&E  RMSS 
Code  B323076A64  S99QAXHB05405  H2590D. 


19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 


Scintillation 
Satellite  Communication 
Propagation  in  Striated  Media 

CO  


ZO'^ABSTRACT  (Continue  on  reverse  side  H necessary  and  identify  by  block  number) 

A series  of  multiple  phase  screen  computations  covered  frequencies  from'-^ 
300  MHz  to  8 GHz,  rms  electron  concentration  fluctuations  from  10  to  10 
electrons/^m^  for  layers  150  km  and  500  km  thick  at  mean  ranges  of  300 
and  1000  km  respectively.  The  Integrated  phase  shift  irregularities  were 
characterized  by  both  exponential  and  inverse  cube  power  spectral  density 
functions  with  various  characteristic  lengths.  When  fading  (scintillation) 
was  significant,  the  Important,  low  amplitude  portion  of  the  receiver 


DD  I 1473  EDITION  OF  I NOV  S5  IS  OBSOLETE 


UNCLASSIFIED 


security  classification  of  this  page  flFhen  Daf«(Enl»r»d) 


UNCLASSIFIED 


security  classification  of  this  PAGEfl«i«n  Dmit  Enlend) 


20.  ABSTRACT  (Continued) 

' signal  is  nearly  Rician.  The  phase  screen  perturl^tions  that  are  associated 
with  structure  of  spatial  wavelengths  less  than  nif  is  the  Fresnel  dis- 
tance) are  those  responsible  for  the  multipath  scintillations.  The  field 
strength  decorrelation  distance  in  the  receiver  plane  is  approximately  i^ times 
the  reciprocal  of  the  rms  transverse  derivative  of  phase  in  the  phase  screens 
when  scintillation  is  significant. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGEf"Tien  Dalm  Enlertd) 


I 


PREFACE 


The  author  would  like  to  thank  Capt.  Leon  A.  Wittwer,  Air  Force 
Weapons  Laboratory,  for  making  available  to  us  his  multiple  phase  screen 
propagation  simulation;  and  Mark  Frolli,  Mission  Research  Corporation,  who 
performed  the  computations  reported  here. 


1 


1 INTRODUCTION 

2 Smil  SCALE  THEORY 

2.1  Environment 

2.2  Dimensional  Scaling 

2.3  Amplitude  Fluctuations 

2.4  Fading  Rates 

3 PHASE  SCREEN  CHARACTERISTICS 

3.1  Relationships  of  Statistical  Descriptors 

3.2  Chesnut  Distribution 

3.3  Natural  Ionosphere 

4 WAVE  PROPAGATION  COMPUTATIONS 

4.1  Multiple  Phase  Screen  Simulation 

4.2  Detrending 

4.3  Grid  Spacing 

4.4  Cases  Computed 

5 RESULTS  OF  COMPUTATIONS 

5.1  Amplitude  Fluctuations —Natural  Ionosphere 

5.2  Amplitude  Fluctuations —Chesnut  Distribution 

5.3  Decorrelation  Distance  — Natural  Ionosphere 

5.4  Decorrelation  Distance — Chesnut  Distribution 

5.5  Decorrelation  Distances —Weak  Scintillation 

5.6  Electric  Field  Magnitude  Autocorrelation  Function 


I 


CONTENTS  (Continued) 


Section 

Page 

6 INTERPRETATION  AND  CONCLUSIONS 

62 

6.1 

Effective  Phase  Variance 

62 

6.2 

Decorrelation  Distances 

64 

6.3 

Shape  of  Autocorrelation  Function 

66 

6.4 

Sumnary 

67 

6.4.1 

Effective  Phase  Shift 

67 

6.4.2 

Decorrelation  Distances 

67 

6.4.3 

Shape  of  Autocorrelation  Function 

67 

6.4.4 

Wavelength  Dependence  of  Effects 

68 

6.4.5 

Research  Requirements 

68 

REFERENCE 

s 

69 

DISTRIBUTION  LIST 


3 


LIST  OF  FIGURES 


Figure 

Page 

1 

Propagation  computation  geometry 

8 

2 

Geometry  of  propagation  problem 

13 

3 

Fade  probability  variation  with  specular-random 

ratio 

19 

4 

Frequency  transmission  function  of  detrend  5-10 

filter 

35 

5 

Field  strength  cumulative  probability  — natural 

ionosphere 

40 

6 

Field  strength  cumulative  probabil ity — natural 

ionosphere 

41 

7 

Effective  phase  variance — natural  ionosphere 

43 

8 

Field  strength  cumulative  probability — Chesnut 

distribution 

44 

9 

Field  strength  cumulative  probabil  ity — Chesnut 

distribution 

46 

10 

Field  strength  cumulative  probability  — Chesnut 

distribution 

47 

11 

Effective  phase  variance  — Chesnut  distribution 

49 

12 

Effective  phase  variance  — Chesnut  distribution 

50 

13 

Decorrelation  distance  — natural  ionosphere  distribution 

51 

14 

Decorrelation  distance  — Chesnut  distribution 

53 

15 

Decorrelation  distance — Chesnut  distribution 

54 

16 

Decorrelation  distance  — Chesnut  distribution 

55 

17 

Decorrelation  distance  as  a function  of  Fresnel 
for  small  phase  fluctuations 

distance 

57 

18 

Autocorrelation  function  of  detrended  field  strength- 
natural  ionosphere 

59 

19 

Autocorrelation  function  of  detrended  field  strength  — 

Chesnut  distribution 

61 

4 


SECTION  1 
INTRODUCTION 


The  detonation  of  nuclear  devices  at  high  altitudes  (E-  and 
i-- layer  altitudes)  produces  copious,  persistent  ionization  which  typically 
divides  into  striations  aligned  with  the  geomagnetic  field.  The  propaga- 
tion of  radio  waves  through  this  medium  results  in  perturbations  of  the 
wave  front  that  may  lead  to  a scintillating,  fading  signal  at  the  receiver. 

A previous  report^  developed  an  analytic  treatment  quantifying 
the  signal  characteristics  if  dimensions  of  the  structure  were  small  com- 
pared to  a Fresnel  zone.  This  work  was  extended  to  treat  fading  rates^ 
under  similar  limitations. 

Unfortunately,  it  appears  that  at  satellite  communication  fre- 
quencies, the  assumption  that  all  of  the  structure  is  small  is  invalid. 
Analytic  techniques  have  not  been  developed  for  structures  about  the  size 
of  a Fresnel  zone.  (Larger  structures  can  be  treated  by  ray  optics  tech- 
niques.) However,  it  is  practical  on  modern,  high-speed  computers  to 
represent  the  striated  medium  as  a series  of  thin  phase  shifting  screens 
and  compute  the  propagation  of  waves  through  this  by  Fourier  propagation 
techniques.  This  technique  is  called  a Multiple  Phase  Screen  (MPS)  compu- 
tation and  produces  a representative  electric  field  distribution  on  the 
receiver  plane. 

A series  of  MPS  computations^  was  run  for  an  L-band  signal  trav- 
ersing an  environment  described  by  an  exponential  power  spectr?  density 
(PSD)  of  integrated  electron  concentration  fluctuations.  The  suits  were 


5 


interpreted  in  terms  of  an  effective  variance  of  phase  shift  integrated 
along  lines  parallel  to  the  line  of  sight  to  the  satellite  transmitter. 

The  effective  variance  was  defined  as  that  which  would  lead  to  a prediction 
of  similar  signal  characteristics  when  used  with  analytic  expressions  devel- 
oped for  small  structure.  The  effective  and  actual  phase  variances  were 
related  through  the  fraction  of  the  PSD  that  corresponded  to  structure  with 
a spatial  wavelength  smaller  than  about  three  times  the  Fresnel  length. 

Methods  of  scaling  the  results  of  the  MPS  computations  to  other 
wavelengths  and  geometries  were  discussed  in  Reference  3.  However,  it 
appeared  desirable  to  obtain  results  over  a wider  range  of  parameters  and 
for  different  phase  shift  power  spectral  densities.  This  report  describes 
the  results  of  MPS  computations  at  frequencies  from  300  MHz  to  8 GHz  for 
two  different  phase  PSD's,  and  at  structural  intensities  such  that  the 
receiver  plane  statistics  varied  from  trivial  perturbations  to  fully  devel- 
oped scintillations.  Primary  emphasis  was  upon  those  environments  for 
which  the  scintillations  would  be  significant  to  a receiver.  Altogether, 
about  2000  simulations  were  run. 
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SECTION  2 

SMALL  SCALE  THEORY 


This  section  will  present  relevant  background  about  the  effect  of 
small  scale  structure  extracted  from  previous  work.  There  will  be  informa- 
tion concerning  the  description  of  the  environment,  dimensional  scaling, 
amplitude  scintillation,  and  fading  rates. 

2.1  ENVIRONMENT 

The  general  configuration  is  that  depicted  in  Figure  1.  An  undis- 
turbed electromagnetic  wave  is  incident  on  a region  of  varying  index  of 
refraction.  The  wave  front  is  assumed  planar,  but  the  extension  to  a spher- 
ical shape  is  only  mathematically  more  tedious  rather  than  essentially  dif- 
ferent. The  wave  propagates  through  the  structured  medium  followed  by  an 
unstructured  one.  The  character  of  the  wave  front  at  a receiver  plane  is 
the  desired  output.  The  assumption  of  a plane  wave  and  infinitely  long 
striations  parallel  to  the  incident  wave  front  renders  the  problem  two- 
dimensional  rather  than  three  and  is  a reasonable  representation  of  the 
nuclear  degradation  problem. 

If  the  striations  are  inclined  to  the  wave  front,  this  analysis 
is  valid  so  long  as  the  inclination  is  not  too  close  to  tt/2  and  if  proper 
care  is  taken  of  the  increased  path  length  through  the  striations. 

The  index  of  refraction  of  a plasma  for  microwaves  (in  regions 
where  the  wave  frequency  greatly  exceeds  the  collision  frequency)  is  given 
by 
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(3) 

(4) 


if  X is  in  centimeters.  For  conditions  of  concern  to  this  study,  the 
index  of  refraction  is  nearly  unity  so  that 


p = 1 - n /2n  . (5) 

e c 

It  is  assumed  that  the  medium  can  be  described  by  a mean  electron 
content  , which  varies  slowly  with  z , and  a fluctuating  increment 

n^(x,z)  such  that 

n.j,  = n^Cz)  + n^Cx.z)  . (6) 


The  phase  advance  compared  to  the  mean  is 


For  validity,  the  thin  phase  screen  approximation  requires  the 
deflection  of  a wave  traversing  a screen  to  be  small  compared  to  structural 
dimensions  of  the  striations.  If  this  is  true,  propagation  through  a phase 
shifting  screen  at  the  middle  of  the  medium  with  a phase  shift  (t>(x)  is  equi- 
valent to  propagation  through  the  actual  medium.  Multiple  phase  screen 
computations  are  applicable  if  the  medium  can  be  divided  into  a series  of 
segments  such  that  each  can  validly  be  considered  thin,  but  the  total  medium 
is  not  thin.  Of  course,  if  the  total  medium  is  also  thin,  an  MPS  computa- 
tion is  still  valid,  but  is  an  unneeded  complication  and  requires  needless 
computing. 


The  phase  shift,  (j),  can  be  any  arbitrary  function.  For  the 

scattering  case  considered  here,  it  is  a random  function  and  it  is  possible 

that  only  its  statistical  description  is  known.  For  this  problem,  (|>  is 

taken  to  be  gaussianly  distributed  because  the  sight  path  intersects  many 

randomly  located  striations.  The  central  limit  theorem  then  leads  to  the 

gaussian  distribution.  The  characteristics  of  the  phase  screens  are 

2 

described  by  the  phase  variance,  , and  either  the  PSD  or  its  equivalent, 
the  autocorrelation  function  R^(5).  The  PSD  and  are  related  by  being  a 
Fourier  transform  pair. 


For  some  purposes  (one  being  the  determination  of  fading  rates), 
it  is  necessary  to  know  the  properties  of  the  transverse  spatial  derivative 
of  the  phase  (|)'.  In  the  appendices  of  Reference  3 it  is  shown  that 


and 


PSD^(k)  dk 


PSD. (k)  dk  . 


(10) 


(11) 
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Equivalently,  these  may  be  obtained  from  the  autocorrelation  function.  If 
we  use  to  denote  the  normalized  autocorrelation  function  of  cf) 

(i.e.,  R^(0)  = 1), 


-o 


:&] 


dc‘-JC=o 


(12) 


)'  is  important  because  the  local  deflection  of  a wave,  0,  is  given  by 


e = ^ 4>' 

2ir  ^ 


(13) 


so  that 


2 _ 2 

*^9  ■ . 2 °4>' 
4Tr  ^ 


(14) 


At  some  point  it  is  necessary  to  relate  the  phase  shift  to  plasma 
parameters.  One  way  that  this  has  been  done  (as  in  Reference  1)  is  to 
visualize  the  medium  as  containing  axisymmetric  striations  with  an  electron 
concentration : 


2 2 

\ " "o  /2a  ) 


(15) 


where  n^  may  be  a function  of  a , r is  the  distance  from  the  striation 
axis,  and  there  is  a distribution  of  striation  sizes  given  by  a probability 
distribution  function,  P(a).  The  following  notation  will  be  used: 


k j 
Uq  a 


k 

nQ(a)  a^  P(a)  da 

0 


(16) 


For  the  normal  case  that  the  striations  overlav  little,  these  expressions 
lead  to  a variance  of  the  local  electron  concentration  of 


2 2 
IT  m n^  a 


(H) 


where  m is  the  number  of  striations  per  unit  area. 
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i 


If  we  define 


and 


a = 
e 


2 

^f  = 


2 2 


2 

no  a 


2 3 
no  a 


2 

no  a 


(18) 


(19) 


the  phase  variance  is 


Jl  _ 

<l>  " .2„2 


X n„  sin  a 
c 


2 3 
no  a 


(20) 


where  a is  the  angle  between  the  mean  direction  of  propagation  and  the 
striation  axes.  In  terms  of  the  variance  of  electron  concentration,  this 
is 


2 2 

- _ 5/2  , a^ 

2 _ 277  L f n_ 

d)  ^ ,2  . a 2 
^ X sin  a e n 


(21) 


The  normalized  autocorrelation  function  is  given  by  the  integral 


%«)  ■ 


/' 

•'-oo 


n?  a^  exp(-^^/4a^)  P(a)  da 


(22) 


2 3 
no  a 


Expression  22  can  be  used  in  Equation  12  to  yield: 

,2 


^5/2  , 

77  L n 

,2  2 

X a sin  a n 
e c 


(23) 
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2.2 


DIMENSIONAL  SCALING 


If  the  thin  screen  approximation  is  valid  (i.e.,  transverse  deflec- 
tion of  energy  is  small  within  the  medium),  the  Fresnel  integral  that  gives 
the  distribution  of  the  field  at  the  receiver  plane  is 


E(y) 


1 

Xr 


/•OO 

-OO 


iTT(x  - 
Xr 


iT 


i<J>(x) 


dx 


(24) 


where  the  geometry  is  shown  in  Figure  2.  Basic  assumptions  are  an  incident 
plane  wave  and  a medium  that  produces  only  moderate  scattering. 


Incident 
Plane  Wave 


Receiver  Plane 


Phase  Screen 


Figure  2.  Geometry  of  propagation  problem. 


Obliquity  factors  of  the  order  of  cos  0 have  been  set  equal  to  unity  so 
that  scattering  at  angles  in  excess  of  0.1  radians  or  so  is  improperly 
treated.  However,  if  scattering  were  that  great,  most  systems  would  not 
operate  satisfactorily.  This  geometry  is  also  applicable  to  propagation 
between  phase  screens  in  the  MPS  approach  where  "receiver  plane"  is  inter- 
preted as  the  next  phase  screen. 
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Inspection  of  Equation  24  indicates  that  there  is  a characteristic 
dimension  associated  with  the  problem:  Ar  . This  is  called  the  Fresnel 
distance: 

= At  . (25) 

The  Fresnel  distance  is  that  distance  along  the  phase  screen  away  from  the 

direct,  straight-line  path  which  results  in  a TT-radian  phase  retardation 

in  the  integral  of  Equation  24.  It  is  basically  a region  of  coherence. 

2 

Outside  that  region,  the  phase  factor  ttx  /Xr  results  in  more  and  more 
closely  spaced  oscillations  of  the  integrand  so  that  the  contribution  to 
the  integral  becomes  small.  If  a new  set  of  variables 


= x/Z^  I 

and 

is  introduced.  Equation  24  can  be  written  as 

/OO  P 

exp  -iTT(x^  - y^)^  + i4»(Xj)1  dx^  , (26) 

-OO  ^ J 

which  is  the  scaled  form. 

The  implication  of  Equation  26  is  that  if  the  phase  shift  at  the 
scaled  location  in  the  phase  screen  x^^  is  maintained  constant  as 
changes,  the  field  at  the  scaled  point  y^  in  the  receiver  plane  remains 
constant.  Thus,  if  the  structure  of  the  phase  variations  is  characterized 
by  a decorrelation  length  Z and  if  and  Z/Z^  are  held  constant,  the 

statistics  in  the  receiver  plane  should  remain  constant  except  that  the 
spatial  correlation  distance  should  vary  as  . 
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2.3 


AMPLITUDE  FLUCTUATIONS 


The  mean  value  of  E can  be  obtained  by  averaging  E over  all 
distributions  of  (J)  . If  the  structure  of  (J)  is  small  compared  to  , 
this  can  be  accomplished  by  multiplying  the  value  of  E as  given  by  Equa- 
tion 2A  by  the  probability  distribution  of  (|)  and  integrating  over  the 
range  of  ({)  . For  the  small  scale  case,  the  order  of  integration  over  x 
and  (})  can  be  inverted  to  yield 

E = Eq  exp (-0^/2)  C27) 

where  E^  is  the  value  of  the  field  in  the  absence  of  any  phase  perturba- 
tions. This  E is  referred  to  as  the  specular  component. 

The  output  of  the  MPS  computations  will  be  interpreted  in  terms 
of  an  effective  phase  variance  crg££  • One  of  the  means  of  estimating  this 
will  be  through  analysis  of  E/E^  where  the  inverse  of  Equation  27  will 
be  used: 

°eff  " ^ '^"(Eq/E)  . (28) 

It  was  argued  in  Reference  1,  that  the  mean  power  is  unaffected 

2 

by  scattering  and  that  power  removed  from  the  specular  component  (l-exp(-a^)) 
is  evenly  divided  between  random  components  in  phase  and  quadrature  to  the 
specular  component.  These  are  also  referred  to  as  real  and  imaginary  com- 
ponents. Thus  the  variances  of  these  are 

= ol  = EQ[l-exp(-o2)]/2  . (29) 


Thus  another  way  of  estimating  an  effective  phase  variance  is 


ohl  = - 2oJ/eJi 


(30) 
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If  there  is  some  large  scale  phase  structure  that  modulates  the 
phase  on  the  receiver  plane  without  affecting  the  amplitude,  it  can  lead  to 
an  apparent  increase  in  the  random  quadrature  component.  Let  the  signal  be 
represented  by 


S = [(E  + g^ap  + ig^a^]  exp(iY)  (31) 

where  g^  and  g^  are  each  random  numbers  from  a unit  variance  gaussian 
distribution  and  y is  the  random  phase  modulation.  Based  upon  the  assump- 
tion that  the  distribution  of  y is  sufficiently  small  that  quartic  terms 
in  Y can  be  dropped  and  that  and  are  equivalent,  the  variance 

of  the  real  component  of  S is 


and  that  of  the  imaginary  component  is 


(32) 


2 2 2 2 

a = a.  + E sin  y . (33) 

S.  1 0 

1 

Thus,  the  variance  of  the  imaginary  component  is  degraded  by  a phase  modula- 
tion to  a greater  extent  that  that  of  the  real  component.  Consequently,  in 
the  data  analysis,  the  variance  of  the  quadrature  component  will  not  be 
used. 


Consider  next  the  amplitude  of  the  signal  |s|  . If  the  expres- 
sion for  S is  taken  from  Equation  31, 

|s|  = ^(E  + g^ap^  + gja^  . (34) 

In  the  weak  scintillation  case,  a.  « E„  , so  that 

1 0 
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2^2 

_ g-O- 

|S|  = E + g a + 

2E 


The  variance  of  IS]  is 
2 


2 1 °i 

"iSl  ' °r  *1? 


(35) 


(36) 


Consequently,  for  cases  of  weak  scintillation,  the  variance  of  the  magnitude 

of  the  field  is  a good  approximation  to  the  variance  of  the  in-phase  random 

2 

component  and  can  be  used  in  Equation  30  to  estimate  '^^ff  ■ 


A parameter  often  used  for  measuring  scintillation  is  the 
index.  Its  square  is  the  variance  of  the  power 


(37) 


For  the  small  structure  case,  E is  Rayleigh-distributed,  and 


or 


= 1 - exp (-20^)  , 

a^ff  = -ljln(l-s2)  . 


(38) 


(39) 


2 

The  ojoblem  with  the  use  of  this  measure  is  that  the  S,  index,  which 

4 

involves  a mean  fourth  power  of  the  signal  amplitude,  is  quite  sensitive  to 

large  values  of  field  strength.  Large  scale  focusing  tends  to  produce  more 

points  of  high  field  strength  than  would  occur  in  a Rayleigh  distribution, 

2 

which  distorts  the  data  somewhat.  In  practice,  values  of  in  excess  of  1 
occurred  for  several  simulations.  This,  of  course,  precludes  the  straight- 
forward application  of  this  measure  to  all  the  data  generated. 


17 


A final  means  of  establishing  an  effective  phase  variance  is  to 
count  the  probability  of  deep  fades.  This  has  the  beauty  of  being  exactly 
the  phenomenon  of  concern.  The  probability  of  very  deep  fades  (greater 
than  20  dB)  is  sufficiently  small  that  the  statistics  on  their  occurrence 
are  not  very  precise  using  the  simulation.  Fortunately,  when  deep  fades 
become  prevalent,  the  probability  density  distribution  is  nearly  uniform, 
or  at  least  has  a constant  gradient  near  the  origin.  For  these  conditions, 
the  probability  of  an  electric  field  amplitude  less  than  a specified  value 
is  proportional  to  the  square  of  that  value.  This  allows  the  extrapolation 
of  the  probability  of  a rare  event  from  a larger  data  base. 


The  specific  means  developed  for  establishing  an  effective  phase 
variance  is  as  follows.  The  probability  of  a 20-dB  or  greater  fade  (ampli- 
tude less  than  0.1)  is  obtained  by  taking  the  probability  of  a fade  to  less 

2 

than  an  amplitude  of  dividing  by  100CE^/Eq)  . Data  is  available 

on  the  probability  of  a 20-dB  fade  as  a function  of  3 (the  specular-to- 
random  ratio)  for  a Rician  distribution  (Figure  3).  Thus,  6 is  obtained 
from  the  extrapolated  probability  of  a 20-dB  or  greater  fade.  3 is 
related  to  the  effective  phase  variance  by 


(40) 


2 

The  problem  with  estimating  this  manner  is  that  it  can 

be  applied  only  over  a limited  range  of  • To  obtain  good  estimates 

with  the  limited  precision  of  fade  probability  requires  that  3 be  between 
about  0.5  and  5,  which  implies  a variation  in  effective  rms  phase  of  2. 
Consequently,  because  steps  of  a factor  of  3 were  used  in  in  most  sets 

of  simulations,  few  good  estimates  were  obtained  for  each  value  of  £ . 
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Specular- to-Random  Ratio 


2.4  FADING  RATES 

Several  potential  causes  of  signal  fading  due  to  propagation 
through  a striated  environment  have  been  investigated:  relative  motion  of 
the  direct  sight  path  and  the  striated  plasma;  relative  motion  of  various 
portions  of  the  plasma;  motion  of  the  source  toward  the  receiver;  the  addi- 
tion of  more  striated  media;  and  decay  of  ionization.  The  primary  cause 
of  fading  here  is  relative  motion  of  the  sight  path  and  the  plasma.  Another 
way  of  visualizing  this  type  of  fading  is  to  consider  it  as  convection  of  a 
nearly  time-invariant  field  past  the  receiver  or  vice  versa.  This  source 
of  fading  is  the  only  one  to  be  considered  here. 

There  are  several  ways  to  estimate  the  decorrelation  distance  of 
the  electric  field  strength  in  the  receiver  plane.  The  variation  in  dis- 
tance to  points  of  constant  phase  in  the  phase  screen  was  used  in  Reference 
4;  the  autocorrelation  function  of  the  electric  field  amplitude  was  devel- 
oped in  Re'^erence  2;  the  beat  distance  for  rays  converging  at  ±a  can  be 
used;  and  the  autocorrelation  function  of  the  electric  field  can  be  used 
in  some  cases.  The  primary  measure  of  fading  rate  or,  equivalently,  fading 
distance  will  be  a decorrelation  distance,  on  the  receiver  plane 
because  this  can  be  obtained  from  the  MRS  computations. 

If  the  striated  environment  is  sufficiently  intense  to  produce 
scintillations,  there  will  be  multi-path  conditions  and  wavelets  inter- 
acting at  the  receiver  plane  which  are  distributed  in  angle  with  a standard 
deviation  . Two  waves  converging  at  angles  of  iOg  will  provide  a beat 
distance 


One  would  expect  the  beat  distance  to  be  somewhat  greater  than  the 
decorrelation  distance  — larger  perhaps  by  a factor  of  two  or  three. 
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The  value  of 
is  related  to 
23; 


"0 

a. 


can  be  obtained  from  Equation  14  and  by  noting  that  0^, 
by  a scale  distance,  s,  obtained  from  Equations  21  and 


= 0^/^  a^ 


(42) 

(43) 


Thus, 


L = a Jo.  . 
b f <p 


(44) 


One 


therefore  expects  the  decorrelation  distance  to  be  a few  times  • 


Consider  the  alternate  derivation  of  the  electric  field  autocorre- 
lation function.  This  is  given  in  Equation  77  of  Reference  1 as: 


R^(C)  = exp{-aj[l  -R^(a]} 


(45) 


Once  the  wave  leaves  the  phase  screen,  the  angular  power  spectral  density 
and,  hence,  the  electric  field  autocorrelation  function  are  independent  of 
the  distance  from  the  screen.  The  character  '’f  the  field  dois  change, 
however.  At  the  screen  the  field  amplitude  is  constant  and  the  decorrela- 
tion is  due  to  phase  modulation.  At  sufficiently  great  distances  and  for 
sufficiently  great  , wave  mixing  occurs,  converting  the  phase  modulation 
to  amplitude  scintillation  with  a statistically  random  phase  or,  equiva- 
lently, into  statistically  independent  in-phase  and  quadrature  components. 

In  this  limit,  the  autocorrelation  function  of  each  of  the  random  components 
becomes  equal  to  that  of  the  initial  field  (Equation  45) . 


for  which 


If  one  defines  the  decorrelation  distance  as  that  value  of 
-1  2 

Rg()!,d)  = e and  if  . only  that  portion  of  R(j,(?) 


where 


21 


it  is  nearly  unity  is  significant.  Consequently,  it  can  be  approximated 
by 


Re(C)  ^1 
2 2 


The  use  of  Equation  45  leads  to 


“d  = ' • 


As  expected,  this  is  somewhat  smaller  than  the  previously  computed  beat 
distance. 


SECTION  3 

PHASE  SCREEN  CHARACTERISTICS 


3.1  RELATIONSHIPS  OF  STATISTICAL  DESCRIPTORS 


In  the  literature  dealing  with  scintillation  produced  by  iono- 
spheric striations  there  is  a variability  in  notation.  Although  the  con- 
ventions used  here  will  not  agree  with  those  used  by  all  workers  in  the 
field,  they  are  as  follows. 


The  normalized  power  spectral  density  of  a zero  mean  random 
function  y in  an  interval  L will  be 


PSD  (k)  = 
y^  ^ 2t\L 


I 


L/2 


-L/2 


y(x)  e dx 


(50) 


and  the  autocorrelation  function  will  be 

R;(^)  =<(y(x)y(x+C)> 
r l/2 

= l/L  i y(x)y(x  + C)dx  . (51) 

■'-L/2 

As  given,  the  power  spectral  density  and  autocorrelation  function  are  not 
normalized.  After  the  autocorrelation  function  is  normalized  to  be  one  at 
C = 0 , 


Ry(C)  = r;(C)/o2  . 


(52) 
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The  power  spectral  density  and  autocorrelation  function  form  a Fourier 
transform  pair.  That  is. 


% ® =/’ 


PSD  (k)  e"^*^^dk  . 

, y 


In  Reference  1 it  was  shown  that  the  variance  of  the  various 
derivatives  of  y could  be  obtained  as  the  value  of  derivatives  of  Ry(?) 
evaluated  at  ^ = 0: 


(h/E)) 


df2n  \ y / 


The  derivative  of  the  autocorrelation  function  can  be  obtained  by  operating 


on  Equation  53.  Thus, 


/C 

C 


k ‘ PSD  (k)  dk  . 

, y 


3.2  CHESNUT  DISTRIBUTION 

A plasma  model  often  considered  was  inferred  by  Chesnut  from 

photographic  data  of  nuclear-detonation  plasma  striations.  As  pointed  out 

by  Chesnut,  this  model  may  be  deficient  in  small  scale  structure  which  was 

* 

obscured  by  photographic  grain  size.  In  any  event,  for  this  : 

R^(C)  = (1  + (57) 

The  i used  here  is  related  to  the  a used  above  by  a = 1//T  . 
a is  used  here  to  conform  to  the  original  Chesnut  analysis  and  the  multiple 
phase  screen  propagation  program  of  AFWL. 
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and 


(58) 

C2/2Z  . 

(59) 

The  corresponding  power  spectral  density  is 


PSD^  = C 4Jl^expC-«^kJl) 


where 


(60) 


(61) 


PSD^  is  taken  to  be  one-sided  (i.e.,  the  portion  associated  with  negative 
k has  been  folded  over  and  added  to  that  associated  with  positive  k ) . 


The  variances  of  derivatives  of  (p 
the  power  series  expansion  of  Rj|j(C) 


are  most  easily  obtained  from 


2 2 2 


(62) 


If  only  striations  smaller  than  a particular  size  are  effective 
at  producing  scintillation,  one  may  reasonably  ask  what  is  the  phase 
variance  associated  with  structure  smaller  than  some  threshold.  This  can 
be  obtained  by  integrating  the  PSD  from  the  threshold  value  k^  to  infinity. 
Thus, 


i-»00 

C4£^  exp(-v^k£)dk 
*^1 

2/2  C Jtexp(-/2k^  £)  . 


(63) 

(64) 
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On  the  basis  of  the  dimensional  scaling  by  the  Fresnel  size,  if 
the  concept  of  a cutoff  is  valid,  should  be  inversely  proportional 

to  the  . Consequently,  for  il  » 2.^  the  effective  phase  variance  should 
decline  rapidly  with  increasing  2-/2,^  , This  reflects  the  lack  of  structure 
at  spatial  wavelengths  much  smaller  than  2,  . Measurements  in  the  natural 
ionosphere  have  not  shown  such  a rapid  structural  cutoff  and  it  is  currently 
felt  by  many  that,  indeed,  film  grain  did  obscure  the  fine  detail  in  nuclear 
test  data.  The  main  reason  for  proceeding  with  this  model  is  that  it  does 
present  a rather  extreme  model  that  may  reasonably  represent  a limiting 
case. 


3.3  NATURAL  IONOSPHERE 


Electron  concentration  profiles  have  been  obtained  by  probes  on 
rockets  and  satellites  traversing  the  ionosphere.  Many  of  the  spectra 
display  the  general  character  of  a one-dimensional  power  spectral  density 


PSD 


2 2 

1 + 


(65) 


— an  "inverse  k square"  relationship.  This  will  be  called  the  natural 
ionospheric  model.  The  line  integral  of  electron  content  or  phase  pertur- 
bation will  have  a power  spectral  density  of  the  form 


<P 


(i.kV)^/2 


(66) 


For  this  model,  it  is  easier  to  work  with  the  power  spectral 
density  than  with  the  autocorrelation  function,  which  is  a Bessel  function. 


The  parameter  2,  in  the  PSD  is  called  the  outer  scale  size. 

There  must  also  be  an  inner  scale  size  which  reduces  the  PSD.  below  that  of 

<P 

Equation  66  in  order  that  the  variance  of  derivatives  of  (J)  remain  finite. 
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There  is  no  good  experimental  data  giving  an  inner  scale  because  it  usually 
appears  to  be  below  the  resolution  (tens  of  meters)  of  experiments  performed 
to  date.  For  demonstration  purposes,  take  an  upper  value  of  k beyond  which 


the  PSD  is  assumed  to  be  zero.  Call  this  k 


m 


The  experimental  data  indi- 


cate that  k^i,»l;  the  upper  wave  number  cutoff  affects  the  normalization 

of  the  PSD^  very  little.  We  can  obtain  by  integrals  of  powers  of  k^  and 

the  PSD,  that 
9 


2 

a. , = Zni2  k H)  • 
<P  ^2  m •' 


(67) 


A specific  representation  of  striations  which  is  convenient  to 

-2 

work  with  and  which  yields  the  k one-dimensional  power  spectrum  is  a 
distribution  of  gaussian  profile  striations  in  which  their  location  is  ran- 
dom and  their  size  distribution  is 


'0 


s < s. 


Ps=<;a/s  Sj<s<S2 
S > St 


(68) 


where  the  inner  scale  s^  is  much  smaller  than  the  outer,  S2  . 
dimensional  electron  concentration  power  spectral  density  is 


The  one- 


1 

PSD  (k)  = ^ 

" /if  (S2-spk2 


/ 1,2  2 ,2  2' 

f -k  s^  -k  $2 
i e - e 


(69) 


At  long  wavelengths  (k  « I/S2) , the  PSD  is  constant 


1,..  '’n^"2""l^ 


PSD  (k)  = 


/if 


(70) 


Between  the  two  scales  (I/S2  « k « l/s^ , the  PSD  exhibits  an  approximately 


inverse  k form 


PSD"(k) 


/ir  (S2  - sp 


at  short  wavelengths  (k»l/Sj^),  the  PSD  cuts  off  exponentially 


a^expC-k^s^) 

PSD  (k)  = — ^ . 

Acs  -s)k^ 


The  phase  variance  is 


oj  = C /F  (Sj  + S2)  . 


The  phase  variance  associated  with  structure  having  a k>k^  can 
be  obtained  by  integrating  Equation  66. 

o?  .C/FlTl (74) 


If  the  outer  scale  is  much  larger  than  the  cutoff  size,  k^X,  » 1 , and 


,2  _ C ✓it 

^ %ilk2  • 


SECTION  4 

WAVE  PROPAGATION  COMPUTATIONS 


In  order  to  investigate  the  detailed  character  of  the  electric 
field  produced  by  plasma  having  structure  with  characteristic  dimensions 
comparable  to  the  Fresnel  length,  computer  simulations  were  run  using  a 
multiple  phase  screen  (MPS)  code®  provided  by  Capt.  Leon  Wittwe^  of  the 
Air  Force  Weapons  Laboratory.  For  completeness,  a short  description  of 
characteristics  of  the  simulation  will  be  included.  This  will  be  followed 
in  this  section  by  a description  of  our  means  for  removing  large  scale 
phase  variations  usually  referred  to  as  detrending,  criteria  used  in  setting 
up  phase  screens,  and  a presentation  of  the  matrix  of  cases  computed. 

4.1  MULTIPLE  PHASE  SCREEN  SIMULATION 

A multiple  phase  screen  propagation  simulation  divides  the  dis- 
turbed environment  into  segments,  collapses  the  phase  shift  in  each  segment 
onto  a single  plane  or  screen,  and  propagates  a wave  through  this  idealized 
representation.  It  is  assumed  that  a segment  of  the  striated  environment 
can  be  chosen  sufficiently  thin  that  there  is  little  deviation  of  energy 
propagating  through  that  segment,  and,  hence,  the  integrated  phase  shift 
experienced  traversing  the  segment  can  be  treated  as  though  occurring  at  a 
single  plane.  The  electric  field  incident  on  the  plane  or  screen  is  multi- 
plied by  the  complex  number  representing  the  phase  shift  to  obtain  the 
electric  field  vector  leaving  the  screen.  Propagation  to  the  next  screen 
is  computed  by  Fourier  analyzing  the  field  into  its  angular  components, 
multiplying  these  by  a propagation  vector  (different  for  each  angular  com- 
ponent), and  transforming  the  propagated  field  back  into  its  spatial  form 
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at  the  next  phase  shifting  screen  or  finally  at  the  receiver  plane.  The 
simulation  utilizes  fast  Fourier  transforms  and  our  computations  used  2048 
points  in  each  phase  screen  and  in  the  angular  distributions. 

There  are  several  ways  one  might  choose  the  distribution  of  phase 
on  the  phase  screen.  One  is  to  randomly  locate  striations  having  the 
structural  properties  and  size  distribution  desired  and  integrate  the  phase 
shift  through  this  explicit  representation  of  striations.  Another  way  is  to 
derive  analytically  a statistical  description  of  the  integrated  phase  and 
from  this  description  generate  directly  a specific  realization  of  the  inte- 
grated phase.  This  latter  procedure  is  the  faster  of  the  two  and  for  the 
cases  to  be  treated,  produces  phase  screens  statistically  equivalent  to 
those  produced  by  the  former  procedure. 

The  phase  screens  generated  have  a phase  whose  amplitude  is 
gaussianly  distributed.  This  is  valid  if  every  ray  through  the  screen 
intersects  several  striations  so  that  the  central  limit  theorem  is  appli- 
cable. To  avoid  the  necessity  of  empty  guard  bands  at  the  edge  of  the 
screen  and  large  angular  Fourier  components  associated  with  a difference 
in  phase  shift  at  the  two  edges  of  the  screen  (the  angular  Fourier  analysis 
assumes  the  field  is  cyclic  in  the  width  of  the  screens),  the  phase  screens 
are  made  continuous  from  one  edge  to  the  other.  The  basic  statistic  of 
the  phase  which  is  used  is  the  autocorrelation  function.  However,  this  is 
Fourier  transformed  to  obtain  the  spatial  power  spectral  density.  The 
amplitude  of  each  Fourier  component  is  multiplied  by  a random  phase  (which 
is  different  for  each  specific  realization)  and  the  phase  screen  is  recon- 
structed by  the  inverse  Fourier  transformation  of  the  randomized  components. 

The  simulation  usually  was  run  twenty  times  for  each  particular 
problem  using  different  phase  screen  realizations.  The  code  generates 
various  averaged  descriptions  of  the  electric  field  at  the  receiver  plane. 
These  include  the  distribution  of  the  electric  field  components  in-phase 
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and  quadrature  to  the  undisturbed  field  and  the  total  field  amplitude.  It 
also  gives  the  two-dimensional  distribution  of  the  electric  vector.  That 
is,  it  gives  a table  of  the  number  of  points  N(n,m)  at  which  ■ 0.05 
<Er^-^+0.05  and  - 0. 05  < + 0.05  . Except  for  a normalizing  fac- 

tor, it  is  the  two-dimensional  probability  density  function  of  E.  It  is 
displayed  as  a plot  with  the  value  of  N(n,m)  printed  at  a coordinate  n,m. 
These  plots  will  be  referred  to  simply  as  "distribution  of  E."  Various 
statistical  descriptors  are  computed  — such  as  the  spatial  autocorrelation 
function  of  the  electric  field. 

4.2  DETRENDING 


Often  the  character  of  scintillation  in  a radio  signal  channel, 
either  measured  through  the  ionosphere  or  computed  with  a simulation,  is 
obscured  by  large,  slowly  varying  phase  shifts.  In  actual  receiver  equip- 
ment, this  slowly  varying  phase  is  removed  by  AFC  circuits  and  phase  track- 
ing loops.  An  extreme  example  of  the  effect  is  the  electric  field  just 
after  the  signal  has  traversed  a single  phase  screen  in  a simulation.  In 
this  case,  the  signal  amplitude  is  unchanged,  but  the  phase  may  vary  widely. 

If  for  the  phase  screen  » 2tt  , the  plot  of  distribution  of  E is  a 
circle  of  constant  amplitude  when  plotted  in  terms  of  E^,  E^^,  and 
E^ = 0 = E^  . In  actuality  there  is,  of  course,  no  scintillation,  although 
the  autocorrelation  function  of  E may  have  a very  short  decorrelation  length. 

Detrending  can  be  approached  from  either  the  time  or  frequency 
domain  of  a receiver  or  equivalently  from  either  the  space  or  spatial  fre- 
quency domain  of  the  electric  field  in  the  receiver  plane  of  a propagation 
simulation.  The  approach  used  here  is  that  considering  the  spatial  variation 
of  the  field  vector.  Detrending  uses  data  from  several  points  to  estimate 
a mean  value  of  the  phase  at  some  other  point.  The  estimated  mean  is  sub- 
tracted from  the  total  phase  and  the  remainder  considered  the  random  com- 
ponent of  interest.  In  this  way,  one  can  attempt  to  isolate  the  large  scale 


phase  fluctuations  and  the  rapid  random  components  due  to  multipath  scin- 
tillations . 


One  must  use  some  care  in  detrending.  If  the  data  points  used  to 
estimate  the  mean  phase  are  too  widely  spread,  there  is  little  detrending 
action  because  the  mean  tends  toward  zero.  On  the  other  hand,  if  the 
points  are  too  limited,  the  estimate  of  the  mean  may  always  be  so  close  to 
the  actual  as  to  eliminate  phase  fluctuations.  With  such  tight  phase  track- 
ing, even  a Rayleigh  distribution  can  be  converted  to  one  which  has  a finite 
real  average  component.  Unfortunately,  it  is  difficult  to  predict  what 
properties  will  lead  to  the  best  reduction  of  the  data. 


The  detrend  design  used  here  operates  with  sampled  data  such  as 
produced  by  the  MPS  calculation.  For  a particular  sample  variable  y^  , 
the  mean  large  scale  component  is  taken  to  be 

^2 

^n  ” ? ^i  ^n+i  (76) 

■^2 

and  the  detrended  value  of  y is 


y 


★ 

n 


^n  - 


b 

n 


(77) 


The  variation  of  the  large  scale  component  can  be  represented  as 

a power  series  in  i . If  the  detrend  filter  is  desired  to  make  v*  inde- 

th  ' " 

pendent  of  the  m power  variation  of  the  background, 

^2 

Ei'"a.  = 0 . (78) 

■^2 


Thus,  in  principle,  it  is  possible  to  make  the  detrend  account  for  all  power 
variations  up  to  2i2;  however,  because  there  is  a random,  noise- 1 ike  com- 
ponent superposed  on  the  slowly  varying  component,  such  a filter  would  be 
noisy. 
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I 


i 
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For  the  current  effort,  the  design  has  been  such  as  to  suppress 
powers  through  the  third.  The  odd  power  components  are  eliminated  if  the 
a^^  sequence  is  an  even  function  (i.e.,  a^  = Thus,  we  are  left  with 

the  restrictions 

^a.  = 1 (79) 

-"2 
^2 

Vi^a.  = 0 . (80) 

“ 1 

~^2 

It  was  also  desired  that  the  numerical  implementation  be  simple  and  that 
the  estimate  of  mean  phase  be  independent  of  the  measured  value  of  the  n 
sample. 


The  simplest  detrend  function  meeting  these  requirements  is  that 


0<|i|<ij 

(81) 

Xi 

11 

This  means 


ilbi  + (i2-  i^)b2  = 0.5 

(82) 

and 

^2 

+ X;  i^b2  = 0 , 

i=l  i=ij^+l 

(83) 

which 

can  be  rewritten  as 

il(bi  - b2)  + i2b2  = 0.5 

(84) 

^1  ^2 

Xi^(b  -b2)  + ° • 

i=l  i=l 

(85) 
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The  expression  for  can  be  similarly  rewritten  as 


b 

n 


i=-i 


n+i 


+ b. 


2 

E y 


1 


i=-i. 


n+i 


(86) 


Parameters  used  in  detrending  to  date  are  given  in  Table  1.  This 
very  simple  detrend  function  was  chosen  for  ease  of  implementation,  rather 
than  for  any  optimization  of  operational  characteristics. 


Table  1.  Detrend  filter  parameters. 


^1 

'2 

b2 

5 

10 

0.140 

-0.020 

10 

0.068333 

-0.0091667 

20 

40 

0.033750 

-0.0043750 

20 

30 

0.045898 

-0.013932 

Figure  4 shows  a typical  filter  transmission  function  and  the 
signal  power  spectral  density  for  a single  sample.  The  signal  and  filter 
responses  were  computed  assuming  the  signal  and  filter  weighting  were  con- 
stant between  sample  times.  The  transmission  function  shown,  T(k),  is  the 
square  of  the  Fourier  transform  of  the  spatial  filter  taking  the  unit 
distance  to  be  the  distance  between  samples.  As  can  be  seen,  most  of  the 
signal  power  is  where  |k|  < tt  . The  filter  shown  is  Detrend  5-10.  Basically, 
the  filter  is  matched  to  the  signal  but  has  both  a sharp  notch  about  k=0 
and  oscillations  superposed  on  it.  For  larger  values  of  i,  and  i , the 
notch  becomes  narrower  (although  the  shape  of  it  still  is  a k relation- 
ship), the  ripple  frequency  becomes  greater  and  the  ripple  amplitude  becomes 
less.  The  relative  amplitude  of  the  ripple  is  about  2(bj-b2)  . If  the 
filter  factors  a^  had  been  tapered  to  avoid  discontinuities  in  the  detrend 
filter,  much  of  the  ripple  could  have  been  suppressed. 


The  field  was  strongly  detrended  for  most  of  the  computations. 

This  was  done  for  a variety  of  reasons.  In  many  equipments,  the  AFC  and 
phase-locked  loops  do  provide  a strong  detrend.  Thus,  the  electric  field 
output  by  the  MPS  computation  is  a realistic  representation  of  the  signal 
in  a receiver. 

In  addition,  the  computed  autocorrelation  function  of  the  detrended 
electric  field  is  a good  representation  of  the  autocorrelation  function  of 
the  field  magnitude  scintillations  and  is  relatively  unaffected  by  any  large 
scale  phase  modulation. 

Finally,  when  the  scintillations  are  weak,  the  properties  of  the 
fluctuations  in  the  detrended  signal  are  dominated  by  the  in-phase  random 
component  so  that  interpretation  of  the  output  is  simplified. 

4.3  GRID  SPACING 

In  order  that  the  MPS  simulation  faithfully  represent  the  desired 
propagation  problem,  the  size  of  the  screens  and  the  spacing  of  points  must 
satisfy  several  criteria.  Five  criteria  were  met  for  the  simulations 
reported. * 

The  width  of  the  screens  was  chosen  to  be  at  least  five  times  the 
phase  shift  decorrelation  distance.  This  is  to  assure  a reasonable  statis- 
tical sampling  of  the  phase  distribution. 

The  grid  spacing  must  be  fine  enough  to  represent  the  small  scale 
structure.  In  all  cases,  structure  below  the  resolution  of  the  grid  con- 
tributed less  than  one  percent  to  effective  phase  variance. 


* 

The  derivation  of  the  quantitative  criteria  are  included  in  Reference  3. 
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It  is  also  necessary  that  the  phase  representation  be  a reasonably 
continuous  function.  The  specific  criterion  used  was  that  the  variation  in 
phase  shift  between  adjacent  points  in  the  grid  represent  less  than  a tenth 
of  a wavelength. 

The  grid  spacing  must  also  be  sufficiently  fine  that  the  maximum 
scattering  angle  be  properly  represented.  Too  great  a grid  spacing  results 
in  aliasing  of  energy  in  wrong  directions.  The  spacing  was  required  to  be 
sufficiently  fine  to  correctly  propagate  a wave  at  three  times  the  standard 
deviation  of  the  scattering  angle.  Thus,  less  than  0.1  percent  of  the 
energy  could  be  aliased. 

Because  the  problem  was  set  up  using  an  assumption  of  a cyclic 
screen  with  a basic  period  of  the  width  of  the  phase  screens,  energy  scat- 
tered out  one  edge  of  the  screen  reappears  at  the  other.  Consequently,  to 
prevent  the  energy  from  a single  point  beating  with  itself  in  a coherent 
manner,  the  rms  transverse  deflection  of  a wave  was  limited  to  less  than 
a tenth  of  the  width  of  the  phase  screens. 

4.4  CASES  COMPUTED 

Computations  were  performed  for  frequencies  from  300  to  8000  MHz 

4 7-3 

for  rms  electron  concentration  variations  from  10  to  10  cm  . The 
striated  media  were  taken  to  be  both  500  km  thick  at  a mean  range  of  1000  km 
and  150  km  thick  at  a mean  range  of  300  km. 

For  the  natural  ionospheric  model,  the  characteristic  length  Z 
was  taken  to  be  0.42,  1.25,  and  3.75  km.  According  to  Wittwer®,  1.25  km 
is  typical  of  strongly  disturbed  natural  equatorial  ionospheres.  These 
parameters  furnished  data  for  values  of  Z/Z^  from  0.38  to  35. 
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For  the  Chesnut  model,  H was  taken  to  be  0.3,  1.0,  and  3.0  km. 

In  addition,  there  was  data  available  from  previous  computations^  for  an  Si 
of  0.1  km.  This  provided  values  of  Si/Si^  from  0.2  to  17. 


Not  all  electron  concentrations  were  run  under  all  conditions 
because  the  small  concentrations  with  high  frequencies  and  long  character- 
istic lengths  result  in  no  effect  and  large  concentrations  at  low  frequen- 
cies and  short  characteristic  lengths  result  in  Rayleigh  distributions.  In 
general,  the  computations  at  a particular  combination  of  wavelength,  charac- 
teristic length  and  mean  range  were  carried  out  over  a range  of  concentra- 
tions such  that  scintillation  ranged  from  trivial  to  fully  developed 
Rayleigh. 


: i' 


SECTION  5 

RESULTS  OF  COMPUTATIONS 


The  results  of  the  MPS  computations  are  divided  in  the  following 
sections  into  amplitude  variations  and  field  strength  structure  fading 
characteristics.  Separate  analysis  of  the  data  for  each  phase  structure 
model  is  presented. 

5.1  AMPLITUDE  FLUCTUATIONS  - NATURAL  IONOSPHERE 

Some  typical  data  are  shown  in  Figures  5 and  6.  The  distributions 
vary  from  one  that  exhibits  very  mild  scintillations  (a  3-dB  fade  occurring 
less  than  two  percent  of  the  time)  to  one  that  is  essentially  Rayleigh 
(greater  than  10-dB  fades  ten  percent  of  the  time) . The  distributions  for 
the  cases  of  lesser  disturbance  can  be  fit  reasonably  well  by  gaussians 
except  in  the  extreme  tails  of  the  distributions.  The  two  straight-line 
approximations  in  Figure  5 are  gaussians  with  variances  obtained  from  the 
spatial  autocorrelation  function  of  the  computed  data. 

The  data  for  which  a.  = 3.7  radians  and  Z = 1.25  km  can  be 

$ 

interpreted  in  terms  of  both  through  the  index  and  o^.  Both 

yield  the  identical  values  of  3 x lo  ^ , which  is  somewhat  fortuitous.  The 
data  for  which  = 2.2  radians  and  Z = 0.42  km  can  be  interpreted 
through  the  S^,  and  the  fading  probability.  These  lead  to  values  of 
a^ff/a^  of  0.028,  0.029,  and  0.029.  The  very  good  agreement  indicates 
that  the  distribution  is  nearly  Rician.  However,  this  environment  is  still 
not  very  severe  from  the  standpoint  of  the  probability  of  deep  fades. 


39 


Percent  Cumulative  Probability 


The  data  for  the  two  environments  given  in  Figure  6 exhibit  signi- 
ficant numbers  of  fades.  The  data  for  which  a,  = 20  radians  and 

<P 

a = 1.25  km  can  be  interpreted  from  the  standpoint  of  S and  the  probabi- 

22  ^ -3  -3 

lity  of  deep  fades  to  yield  values  of  1.4  xio  and  0.85x10  , 

respectively.  That  the  index  leads  to  a larger  value  is  typical  because 

of  the  emphasis  in  computing  of  large  amplitude  points. 


2 2 

Figure  7 is  a summary  of  all  the  data  about 


natural  ionospheric  model.  The  points  about  a factor  of  two  above  the 

curve  fit  to  the  data  are  values  deduced  from  the  S.  index.  The  curve  fit 

4 

to  the  data  is 


(87) 


5.2  AMPLITUDE  FLUCTUATIONS — CHESNUT  DISTRIBUTION 


The  next  set  of  figures  illustrates  data  generated  using  the 
Chesnut  distribution.  Figure  8 shows  the  reproducibility  of  the  data.  The 
two  curves  are  each  the  average  statistics  of  20  separate  realizations  of 
the  phase  screens.  The  general  character  of  the  two  curves  is  the  same. 

At  the  1 percent  (or  99  percent)  probability  level,  the  number  of 
smaller  amplitude  points  is  20  per  simulation,  or  400  per  set  of  20  simula- 
tions. Thus,  if  each  point  were  independent,  one  would  expect  a standard 

4 

deviation  of  20  out  of  4 x lo  points,  or  0.05  percent.  At  the  99  percent 
cumulative  probability  point,  the  separation  of  the  two  curves  is  about 
±0.1  percent,  which  indicates  not  great  correlation.  On  the  other  hand,  at 
the  1 percent  probability  point,  the  curves  separate  a large  amount.  Evi- 
dently, the  low  amplitude  points  occur  in  groups. 
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Percent  Cumulative  Probability 


The  other  significant  characteristic  of  the  curves  is  that  they 
are  decidedly  not  Rician.  The  probability  of  low  amplitude  points  is  about 
the  same  as  for  the  data  shown  in  Figure  5 for  the  natural  ionospheric  model 
with  = 3.7  , but  the  probability  of  electric  fields  more  than  1.4  times 
the  mean  value  is  ten  times  greater  — this  is  characteristic  of  data  gener- 
ated for  the  Chesnut  model  when  Z is  considerably  greater  than  8.^  . 
Segments  of  the  phase  screen  large  compared  to  act  as  coherent  lenses. 

In  contrast,  for  the  natural  ionospheric  model,  there  is  sufficient  small 
scale  structure  to  break  up  the  coherency  from  extended  regions. 

The  development  of  the  focusing  is  also  shown  in  Figure  9.  Here 
the  increase  in  rms  phase  changes  the  low  amplitude  portion  of  the  distri- 
bution function  little,  but  the  large  amplitude  portion  markedly.  The  focal 
length  associated  with  the  rms  value  of  the  second  derivative  of  phase  shift 
is  1500  km  for  the  screen  for  which  = 40  radians.  Inasmuch  as  the  mean 
distance  to  the  phase  screens  was  1000  km,  it  is  not  surprising  that  focus- 
ing was  evident.  The  data  making  up  this  curve  also  had  the  property  of  an 

S.  index  in  excess  of  1 and  a deep  fading  probability  corresponding  to  a 

^ 2 2 
0^  of  less  than  0.3  . Thus,  for  these  curves,  depended  to  a 

greater  extent  on  the  type  of  data  being  interpreted. 

Figure  10  shows  some  data  for  a smaller  value  of  Z/Z^  . For  this 
data,  the  low  amplitude  portion  of  the  probability  distribution  is  more 
nearly  Rayleigh.  For  most  systems  we  have  studied  so  far,  the  most  signifi- 
cant effect  of  scintillations  is  loss  of  signal  during  fades.  For  condi- 
tions such  that  fades  are  prevalent,  the  cumulative  probability  distribution 
function  is  nearly  Rayleigh  in  the  all-important  low  amplitude  portions  of 
the  distribution.  This  is  probably  not  too  surprising  because  all  it  implies 
is  that  the  two-dimensional  probability  density  function  in  the  complex 
plane  is  nearly  constant  in  a region  near  the  origin. 
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Percent  Cumulative  Probability 

Figure  9.  Field  strength  cumulative  probability  — Chesnut  distribution 
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The  results  of  the  MPS  computations  were  interpreted  in  terms  of 

2 2 

Oeff/o^  and  the  data  are  shown  in  Figures  11  and  12.  Figure  11  includes 
all  of  the  interpreted  data;  Figure  12  includes  only  those  interpretations 
based  upon  the  probability  of  deep  fades.  It  can  be  seen  that  up  to  a 
value  of  about  3 for  , the  data  can  be  represented  by 

a^ff/a^  = exp(-2.8  £/£f)  . (88} 

For  larger  values  of  , the  data  exceed  the  exponential.  A fourth 

power  law  is  fit  to  the  data  in  Figure  12.  In  Reference  3 it  was  shown 

that  if  the  cause  of  the  variation  in  signal  strength  is  large  scale 

-4 

focusing  rather  than  multipath,  the  (Jl/£^}  relationship  is  appropriate. 
Thus,  for  the  Chesnut  distribution  and  large  values  of  , the  observed 

trend  in  the  data  seems  reasonable. 


5.3  DECORRELATION  DISTANCE — NATURAL  IONOSPHERE 

The  decorrelation  distance  was  taken  to  be  that  distance  along 
the  receiver  plane  such  that  the  autocorrelation  function  of  the  random, 
in-phase  component  of  the  electric  field  was  less  than  the  value  at  the 
origin  by  the  factor  e ^ . In  those  cases  for  which  data  with  different 
detrend  functions  were  available,  that  function  giving  the  tightest  phase 
track  was  used.  Thus,  the  autocorrelation  function  used  was  very  close  to 
that  for  the  electric  field  magnitude.  For  most  of  the  cases,  the  variance 
of  the  quadrature  component  was  much  smaller  than  that  of  the  in-phase  com- 
ponent, after  detrending. 

The  decorrelation  distances  for  simulations  using  the  natural 

ionospheric  model  and  relatively  large  phase  shifts  are  shown  in  Figure  13. 

Although  there  is  considerable  scatter  of  the  data  points,  they  do  cluster 

about  lines  with  the  relationship  > except  at  the  smaller  values 

of  o.  . (The  decorrelation  distance  for  small  values  of  will  be 

<P  (}> 
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Figure  11.  Effective  phase  variance  — Chesnut  distribution 
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Figure  12.  Effective  phase  variance  — Chesnut  distribution 
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Figure  13.  Decorrelation  distance — natural  ionosphere  distribution 
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discussed  more  in  Section  5.5.)  However,  the  decorrelation  distance  at 

small  values  of  a,  is  rather  academic  because  scintillation  is  minor  in 
<P 

those  cases.  Thus,  for  cases  of  significant  scintillation,  the  decorrela- 
tion time  is  the  time  required  for  relative  motion  of  the  plasma  and  sight 
path  to  convect  the  field  distribution  a distance  in  the  receiver  plane 
that  is  inversely  proportional  to  the  rms  phase  shift  of  the  screens. 


2 

If  one  uses  the  data  for  which  the  S.  index  exceeds  0.6,  the  value 

4 

of  is  0.58  for  an  £ of  420  m;  0.54  for  an  £ of  1250  m;  and 

0.55  for  an  £ of  3750  m.  It  should  be  noted  that  for  these  cases  of  sig- 

nificant scintillation,  the  data  varied  about  the  simple  inverse  power  fit 
by  about  ±50  percent  for  all  computed  cases  for  wavelengths  from  3.75  to 
100  cm,  for  outer  scale  lengths  from  420  to  3750  m,  and  for  mean  ranges  to 

the  phase  screen  of  300  and  1000  km.  Thus,  the  data  cannot  be  taken  to  be 

a particularly  pathological  case,  but  must  be  assumed  to  have  a reasonably 
general  applicability. 

5.4  DECORRELATION  DISTANCE  — CHESNUT  DISTRIBUTION 


Similar  decorrelation  distance  data  is  shown  in  Figures  14,  15, 
and  16  for  data  generated  with  the  Chesnut  model  of  phase  screen  phase 


shift.  As  can  be  seen,  for  large  values  of 
the  decorrelation  distance  falls  off  as  a 


-1 

<P 


<t> 


(significant  scintillation). 


The  lines  at  small  values  of  are  merely  to  identify  points 

associated  with  nearly  equal  Fresnel  distances  (given  in  meters),  which  is 
a subject  of  discussion  in  Section  5.5. 


The  particular  lines  drawn  for  approximating  the  data  for  large 
scintillation  cases  are  not  best  fits,  but  are  the  theoretical  lines  given 
by  Equation  48,  as  will  be  discussed  in  Section  6.2.  The  best  fitting  lines 
with  a slope  of  -1  are  equal  to  0.9  when  £ = 100  m;  equal  to  1.32 
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when  a = 300  m;  equal  to  1.59  when  £ = 1000  m;  and  equal  to  1.38  when 
£ = 3000  m. 

5.5  DECORRELATION  DISTANCES — WEAK  SCINTILLATION 

Figure  17  is  a different  plot  of  the  decorrelation  distance  for 
2 

weak  scintillations  (S^  < 0.1).  Here  the  decorrelation  distance  is  plotted 
in  terms  of  the  Fresnel  distance.  It  can  be  seen  in  Figures  14,  15,  and  16 
that  at  small  values  of  , £^  becomes  relatively  constant  for  a fixed 
Fresnel  length. 

All  the  data  for  small  values  of  o^  for  the  natural  ionosphere 

4> 

fit  a single  curve  of  the  form  £^  = 0.324  £^  and  the  variance  of  individual 
computations  (42  cases)  is  only  0.018  £^.  Thus,  for  cases  of  small  scin- 
tillation, the  characteristic  fading  rate  for  the  natural  ionospheric  model 
is  associated  with  a spatial  autocorrelation  distance  of  a thira  of  a 
Fresnel  distance. 

The  decorrelation  for  the  Chesnut  model  is  not  so  uniform.  For  a 
characteristic  length  of  300  m for  the  distribution  of  phase  shift,  the  data 
for  all  wavelengths  and  for  both  300-  and  1000-km  mean  separation  of  the 
phase  shifting  screens  and  receiver  plane  lie  on  a single  curve  with  a slope 
of  about  0.5.  Moreover,  at  Fresnel  distances  in  excess  of  300  m,  the 
decorrelation  distances  are  nearly  the  same  as  for  the  natural  ionosphere. 

One  is  inclined  to  speculate  (without  definitive  data)  that  if  there  is  con- 
siderable structure  in  the  phase  screens  with  chaiacterisi’ ic  dimensions 
smaller  than  the  Fresnel  length,  the  electric  field  magnitude  in  the  receiver 
plane  will  show  a decorrelation  length  comparable  to  the  Fresnel  length. 

The  data  for  the  Chesnut  distribution  and  £ = 1000  m displays  a 
different  character.  It  splits  into  two  branches  at  small  values  of  £^ 
dependent  upon  the  range  to  the  phase  shifting  screens.  For  the  screens  at 
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a mean  range  of  300  km,  the  variation  of  is  nearly  proportional  to  the 

square  of  !L^  . This  implies  that  it  is  directly  proportional  to  the  wave- 
length of  the  propagating  signal.  The  cause  of  this  has  not  been  identified. 
It  should  be  reiterated  that  for  this  distribution  of  phase  shifts  in  the 
screens  and  the  large  ratio  of  characteristic  of  the  data,  the  scin- 

tillation parameter  is  small  and  the  major  effect  is  large-scale  focusing 
rather  than  multipath.  Thus,  one  might  expect  the  decorrelation  distance  to 
exceed  that  associated  with  true  multipath.  In  Reference  3,  it  was  esti- 
mated that  multipath  would  dominate  effects  for  the  natural  ionospheric 
model  and  for  £ = 300  m for  the  Chesnut  model,  but  that  large-scale 
focusing  would  dominate  for  £ = 1000  m and  < 10  (which  is  the  region 
for  which  £^  is  independent  of  • 

In  conclusion,  one  has  the  feeling  that  if  small-scale  fluctua- 
tions in  the  phase  shift  lead  to  multipath  conditions,  the  detrended  signal 
magnitude  will  exhibit  a decorrelation  length  comparable  to  the  FresneJ 
distance  (the  current  simulations  show  £j  = £^/3  ) . 

5.6  ELECTRIC  FIELD  MAGNITUDE  AUTOCORRELATION  FUNCTION 


: The  decorrelation  distance  discussed  in  the  previous  two  sections 

I was  the  distance  for  the  autocorrelation  function  of  the  real  part  of  the 

t detrended  electric  field  to  fall  a factor  of  e ^ . For  the  evaluation  of 

the  effect  of  propagation  through  a gtriated  medium  upon  the  performance  of 
a communication  system,  it  is  often  necessary  to  know  the  shape  of  the  auto- 
correlation function  as  well  as  the  "decorrelation  distance." 

Some  autocorrelation  functions  are  shown  in  Figure  18  for  the 
I natural  ionosphere  model.  Although  there  is  considerable  variability  in 

I the  curves  when  the  autocorrelation  function  is  much  less  than  the  value  at 

the  origin  (i.e.,  the  normalized  autocorrelation  function  is  much  less  than 
unity),  this  portion  of  the  autocorrelation  function  has  very  minimal 
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Autocorrelation  Function 


effects  upon  the  characteristics  of  a sequence  of  values  of  the  function. 
Where  the  normalized  autocorrelation  function  exceeds  perhaps  10  to  20 
percent,  the  autocorrelation  function  is  nearly  quadratic  with  zero  slope 
at  the  origin. 

Autocorrelation  functions  of  the  field  strength  produced  by  the 
Chesnut  model  (Figure  19)  appear  quite  different  unless  the  intensity  of 
the  striations  is  great.  In  this  context,  "great"  implies  that  the  effec- 
tive rms  focal  length  of  the  phase  screen  is  considerably  shorter  than  the 
separation  of  the  phase  screens  and  the  receiver  plane.  This  means  that 
there  is  much  ray  crossing  between  the  phase  screens  and  the  receiver.  If 
the  focal  length  is  comparable  to  or  greater  than  the  separation,  which  is 
the  case  for  the  curves  shown  in  the  figure,  the  autocorrelation  function 
of  the  in-phase  component  of  the  random  signal  tends  to  be  more  nearly 
exponential  with  some  indication  of  potential  rounding  at  the  origin. 
However,  the  points  shown  in  the  figure  correspond  to  the  grid  spacing  in 
the  computation,  so  that  there  is  no  additional  information  about  the 
detailed  shape  of  the  curves  close  to  the  origin. 
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SECTION  6 

INTERPRETATION  AND  CONCLUSIONS 


I 

i 

The  data  base  included  in  the  computations  reported  here  is, 
admittedly,  limited.  One  distribution  of  random  phase  shift  has  little  j 

high  frequency  content,  the  other  has  considerable;  but  although  they  do  j 

represent  quite  different  distributions,  there  were  only  two  distributions  | 

used.  This  section  will  interpret  the  results  in  terms  of  the  theory  pre- 
sented earlier  and  some  conclusions  will  be  stated  which,  it  is  hoped,  are  1 

more  generally  applicable  than  to  the  two  distributions  studied. 

! 

6.1  EFFECTIVE  PHASE  VARIANCE  | 

I 

1 

The  effective  phase  variance  for  the  natural  ionospheric  model 
(Figure  7)  was  well  fit  by  the  expression  (Equation  87): 


This  is  exactly  the  same  form  as  that  derived  for  that  portion  of  the  phase 
variance  of  the  phase  screen  that  corresponds  to  structure  of  wave  number 
greater  than  some  cutoff  value  (Equation  74) . Let  us  say  that  at  least  for 
this  distribution  of  phase  shifts,  the  portion  which  is  effective  at  pro- 

imit.  Call  the  wave 
Xj  (=  2Tr/kj)  . If 
and  Equations  74 


= 2/£^ 


(89) 


or 

. (90) 

If  similar  logic  is  applied  to  the  Chesnut  model  data  and  Equa- 
tion 88  (the  experimental  results)  is  compared  to  Equation  64,  the  result 
is  identical  to  that  for  the  natural  ionospheric  model,  • 

Thus,  in  cases  of  significant  scintillation  it  seems  reasonable 
to  compute  an  effective  phase  variance  associated  with  structure  having  a 
spatial  wavelength  shorter  than  , and  to  use  this  effective  phase 

variance  in  the  manner  described  for  small  scale  structure  in  Reference  1. 
That  is,  the  specular  power  is  reduced  over  the  free  field  power  by 

P/Pq  = exp(-a^^^)  (91) 

and  the  remaining  power  is  divided  between  in-phase  and  quadrature  random 
components.  Although  the  division  is  uneven  for  structure  comparable  in 
size  to  if  scintillation  is  not  too  great,  for  cases  of  interesting 

degradation  with  some  detrending  device  (AFC,  for  instance),  little  error 
should  be  introduced  into  simulations  by  taking  the  variances  of  the  two 
components  as  equal,* 

This  result  with  reference  to  effective  phase  variance  also 
implies  that  if  the  propagation  medium  extends  over  a long  path  length, 
then  in  the  integration  of  the  environmental  parameters  along  the  sight 
path  one  should  compute  the  incremental  effective  phase  variance  for 


Some  systems  may  be  more  sensitive  to  the  division  between  random 
signal  components  than  others.  The  actual  numerical  results  from  MPS 
calculations  can  be  employed  directly  in  receiver  simulations  to  assess 
the  sensitivity  to  variations  in  signal  statistical  properties. 
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segments  of  the  path  and  sum  them.  Thus,  the  variation  with  distance  of 
Fresnel  distance  and  structural  dimensions  can  be  accounted  for. 

6.2  DECORRELATION  DISTANCES 


It  is  interesting  to  compare  the  measured  decorrelation  distances 
with  the  theoretical  values  far  from  the  phase  screen  as  given  in  Equation 
48  . For  the  natural  ionospheric  model.  Equation  67  relates 


a,  and  o,, 
<P  <t) 


9 o 

0^,  = -I  Jln(2k  . 


(92) 


Inserting  the  value  of  the  inner  scale  wavelength  in  terms  of  the 

inner  scale  wave  number  yields 


^(4^  j 


1/2 


(93) 


The  appropriate  value  to  use  for  is  unclear.  The  descrip- 

tion of  the  MPS  simulation  (Reference  5)  implies  that  it  is  about  20  m. 

If  it  were  shorter  than  that,  the  highest  frequency  structure  would  not  be 
represented  in  the  simulation  because  most  cases  used  a grid  spacing  of 
10  m.  If  one  requires  two  points  per  cycle  (necessary  to  avoid  aliasing), 

A once  again  comes  out  to  be  about  20  m.  If  this  minimum  wavelength  were 
m 

in  error  by  a factor  of  two,  it  would  make  a difference  of  only  6 percent 

in  the  theoretical  value  of  . Table  2 shows  the  experimental 

value  for  the  mean  value  of  using  data  for  which  > 0.6  and 

also  shows  the  theoretical  values  using  A = 20  m . 

m 


The  theoretical  and  experimental  values  are  sufficiently  close  that  one 
feels  reasonable  confidence  in  computing  in  this  manner.  It  is  inter- 

esting that  from  the  standpoint  of  the  decorrelation  distance,  it  is  the 
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total  value  of  that  is  important,  not  the  effective  value.  The  reason, 

4> 

of  course,  is  the  use  of  only  that  data  associated  with  large  amounts  of 

scintillation.  The  effective  a determines  when  scintillations  start, 

eft 

but  once  they  have  started,  multipath  occurs  for  rays  converging  according 
to  the  total  . Actually  the  straight-line  fits  with  unit  slope  in 

Figure  13  are  the  theoretical  curves  and  can  be  seen  to  represent  the  data 
well  within  the  scatter  of  individual  points. 


Table  2.  Normalized  decorrelation  distances  for  cases  of 
significant  scintil lation —natural  ionosphere. 


£ 

(meters) 

Experimental 

Theoretical 

420 

0.575  ±0.038 

0.60 

1250 

0.543  ±0.028 

0.55 

3750 

0.55  ±0.048 

0.51 

What  is  the  case  for  the  Chesnut  distribution?  In  that  set  of 
computations  is  just  (Equation  62)  . Thus,  theoretically 

/2  . (94) 

This  was  used  to  determine  the  fit  to  the  data  in  Figures  14,  15,  and  16. 

The  fit  is  within  10  percent,  except  for  the  data  for  which  H = 100  , where 
experimentally  is  0.9  rather  than  1.4  . 

It  certainly  appears  reasonable  to  approximate  the  decorrelation 
distance  by  unless  additional  research  indicates  otherwise. 
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6.3 


SHAPE  OF  AUTOCORRELATION  FUNCTION 


For  the  phase  screens  that  appear  to  give  true  multipath  in  the 
sense  of  many  interfering  rays  (the  natural  ionosphere),  the  spatial  auto- 
correlation function  has  zero  slope  and  is  concave  downward  at  the  origin. 
This  leads  to  well  behaved  spatial  derivatives  of  the  electric  field  (and 
well  behaved  time  derivatives  as  it  convects  past  a receiver).  As  the  time 
intervals  between  samples  of  the  signal  decrease,  the  time  derivative  of 
the  signal  becomes  constant,  so  that  the  change  in  the  signal  becomes  pro- 
portional to  the  time  interval. 

In  contrast,  for  some  cases,  the  Chesnut  model  leads  to  nearly 
exponential  autocorrelation  functions.  This  can  be  explained  in  terms  of 
ray  optics.  The  large  scale  focusing  associated  with  the  lack  of  high 
spatial  frequency  components  in  the  phase  screen  results  in  optical  caustic 
curves.  Across  such  a caustic,  the  electric  field  may  change  rather 
abruptly.  It  is  such  discontinuities  in  the  field  that  can  lead  to  an  auto- 
correlation function  which  has  a non-zero  slope  and  an  indeterminate  curva- 
ture at  the  origin. 

In  actual  practice,  diffraction  should  limit  the  abruptness  of 
changes  across  the  caustics  and  hence  round  off  the  autocorrelation  function 
at  the  origin.  Even  in  the  worst  set  of  data  shown  (the  data  in  Figure  19 
for  which  % = 1000  m,  X = 10  cm,  = 58  rad,  and  F = 300  km),  there  is  a 
hint  of  this  round  off.  There  is  one  more  effect  that  rounds  off  the  auto- 
correlation function.  The  signal  actually  detected  by  the  receiver  has  an 
autocorrelation  function  which  is  the  convolution  of  the  field  strength 
autocorrelation  function  and  the  antenna  aperture. 

It  is  the  feeling  of  the  author  that  it  is  unlikely  that  an  actual 
distribution  of  phase  shifts  will  be  so  deficient  in  high  spatial  frequency 
components  as  the  Chesnut  distribution  and  in  fact  that  there  will  always 
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be  a region  near  the  origin  where  the  autocorrelation  function  can  be  fit 
by  a concave- downward  quadratic. 


6.4  SUMMARY 

6.4.1  Effective  Phase  Shift 

The  environmental  phase  shift  effective  in  producing  scintillation 
in  microwave  satellite  communication  systems  can  be  obtained  from  the  PSD 
of  the  integrated  phase  fluctuation 

This  phase  shift  variance  can  then  be  used  as  though  it  came  from  structure 
small  compared  to  the  Fresnel  zone. 

6.4.2  Decorrelation  Distances 

The  decorrelation  distance  of  the  electric  field  in  the  receiver 

plane  is 


I 


d 


C96) 


if  scintillations  are  significant.  ( is  the  variance  of  the  transverse 
derivative  of  integrated  phase  in  the  environment.) 


6.4.3  Shape  of  Autocorrelation  Function 

For  the  natural  ionospheric  model  and  for  many  conditions  for  the 
Chesnut  model,  the  autocorrelation  function  is  well  represented  by  a down- 
ward open  parabola  with  a width  commensurate  with  the  decorrelation  dis- 
tance. In  some  cases,  for  the  Chesnut  model,  large  scale  focusing  leads  to 
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caustic  curves  and  an  autocorrelation  function  which  is  parabolic  over  a 
limited  region  and  then  becomes  exponential. 


6.4.4  Wavelength  Dependence  of  Effects 


It  is  interesting  to  compare  the  value  of  the  variance  of  the 
local  electron  concentration  fluctuations  necessary  to  produce  a given 
level  of  signal  scintillation  (i.e.,  the  onset  of  major  scintillations). 

If  i > , the  effective  phase  variance  is  proportional  to  , where 

n is  about  2 for  the  natural  model  and  4 for  the  Chesnut  model.  Thus,  for 


« X 


1/2 


, so  that  the 


a particular  value  of  ; also, 

2 -n/2  4)  f ^ 2 2 2 

needed  a , « X .On  the  other  hand,  for  a fixed  environment  0,  “ X a 

The  final  implication  is,  then,  that  to  produce  a given  level  of  scin- 
tillation requires  a local  electron  concentration  variance  that  depends  upon 
X for  the  natural  ionospheric  model  and  X for  the  Chesnut.  However, 
at  the  onset  of  scintillations,  the  fading  rate  will  be  faster  at  the 
shorter  wavelengths  because  the  total  will  be  greater. 


6.4.5 


Research  Requirements 


In  order  to  evaluate  satellite  communication  systems  with  pre- 
cision, it  is  necessary  to  develop  the  capability  of  performing  the  inte- 
gral of  Equation  95.  Thus,  the  PSD  of  the  integrated  environmental  phase 
shift  must  be  available  at  least  at  spatial  wavelengths  less  than  typical 
Fresnel  zones.  Current  weapons  effects  phenomenology  does  not  provide  this 
information. 

Another  reason  for  needing  the  research  to  determine  the  PSD  is 
to  allow  a more  definitive  choice  between  the  parabolic  and  exponential 
spatial  autocorrelation  functions  in  the  receiver  plane.  (Although  even 
for  the  Chesnut  distribution,  only  a limited  range  of  problem  parameters 
leads  to  exponential  autocorrelations.) 
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